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(54) Adaptive error diffusion method 

(57) An error diffusion method suitable for produc- 
ing an output image from an input image having a set of 
digitized continuous-tone pixels is disclosed. The 
method includes the st^DS of computing an image activ- 
ity signal; computing a set of activity weights from the 
image activity signal; computing a tittered input value for 
a digitized continuous-tone input pixel responsive to the 
activity weights and computing a filtered output value 
responsive to the activity weights for each of the possi- 
ble output levels. The method further includes selecting 



the output level in response to the filtered input value 
and the filtered output value for each of the possible out- 
put levels according to an error criterion; determining an 
error signal between the filtered input value and the fil- 
tered output value for the selected output level; and 
weighting the error signal and adjusting the filtered input 
values for nearby pixels which have not yet been proc- 
essed. 
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Description 

Cross Reference to Related Appiication 

Reference is made to commonly assigned U.S. 
Serial No. 08/399,678 filed 3/7/95 to Kevin E. Spaulding 
and Douglas W. Ckjuwenhoven, entitled IMPROVED 
ERROR DIFFUSION METHOD, the teachings of which 
are incorporated herein. 

Held of the Invention 

The present invention relates to digital image 
processing and more particularly to a method tor digit- 
ally halftoning a continuous-tone image using error dif- 
fi^ion. 

Background of the Invention 

Digital halftoning is a technique employing digital 
image processing to produce a halftone output image 
from a continuous-tone input image. In the digital half- 
toning technique, a continuous-tone image is sampled, 
for exarrple, with a scanner and the samples are digi- 
tized and stored in a computer. The digitized samples 
(or "pixels**) consist of discrete values typically ranging 
from 0 to 255- To r^rcduce this image on an output 
device capable of printing dots of one gray level (e.g. 
bl . jk) it is necessary to create the sensation of multiple 
gray levels by suitably distributing the printed dots in the 
output image. The is accomplished by converting the 
continuous-tone image to a binary output image using 
some form of halftoning technique. 

One prior art method of digital halftoning is kr>own 
as error diffusion. FIG. 1 shows a block diagram 
describing a basic error diffusion technique. The contin- 
uous-tone input value for column i. and row j of the input 
image is given by yjj. For purposes of Illustration it will 
be assumed that the continuous-tone input values span 
the range from 0 to 255. The continuous-tone input 
value for the current input pixel is thresholded 1 0 to form 
the output value bj j. The threshold operator will return a 
0 for any continuous-tone input value below the thresh- 
old, and a 255 for any continuous-tone input value 
above the threshold. A difference signal generator 12 
receives the continuous-tone input value and the output 
value, and produces a difference signal representing 
the en'or introduced by the thresholding process. The 
difference signal is multiplied by a series of error 
weights using a weighted error generator 1 4, and is pro- 
vided to an adder 16 which adds the weighted differ- 
ence signal to the continuous-tone input values of 
neaiby pixels which have yet to be processed to form 
modified continuous-tone input values. The propagation 
of the errors made during the quantization process to 
the nearby pixels insures that the arrtiimetic mean of the 
pixel values is presen/ed over a local image region. 



FIG. 2 illustrates a typical set of error weights 14 
which can be used to distribute the errors to the nearby 
pixels. 

FIG. 3 shows a typical image generated using this 

5 simple error diffusion technique. 

An artifact that is typically associated with error dif- 
fusion halftoning techniques is known as Vorms,"* 
Worms are formed when the black or white output pixels 
appear to string together in an area which should be 

70 othenftrise uniform. Worm artifacts can be clearly seen 
in several areas of the sample image shown in FIG. 3. 
such as the light and dark ends of the gray wedge. Sev- 
eral of these worm artifacts are labded as 1 8A-C. Many 
modifications to the basic error diffusion technique have 

75 been described which attempt to eliminate these worm 
artifacts. One such modification has been described by 
Goertzel et al in U.S. Patent No. 4.654.721 issued 
March 31, 1987. In their disclosure, each input pixel is 
represented by a block of , output pixels. The appearance 

20 of worms is avoided by employing a rarKiom number 
generator to determine the distribution of error between 
two adjacent blocks of pixels, such that a random frac- 
tion of error is distributed to one block and the remain- 
der of the error is distributed to the other. In another 

25 technique, sometimes called error diffusion with dither, 
the error criterion is randomly varied to reduce the 
appearance of worms. Altinough these techniques, and 
others like them, reduce the appearance of worms, they 
also ino-ease the random noise in the output image. To 

30 the extent that this random noise occurs in the frequen- 
cies visible to the human viewer, the image is tiiereby 
degraded. 

A more desirable approach has been disclosed by 
Sullivan (U. S. Patent 5.051,844 issued September 24. 

35 1 991). This method will be referred to as visual error dif- 
fusion. Sullivan teaches the use of a visual filter derived 
from the response of the human visual system to com- 
pute a visually perceived output value. The output level 
is chosen which gives the smallest error between the 

40 continuous-tone input value and the visually perceived 
output value. The error signal which is propagated to 
the nearby image pixels is then computed by taking the 
difference between the continuous-tone irpul value and 
the visually perceived output value rather than the out- 

45 put value itself. FIG. 4 illustrates a flow diagram for this 
method. The simple threshold 10 in the conventional 
error diffusion method has been replaced by a selector 
30 which determines the output pixel value bj j. The 
selection by the selector 30 is made by using a causal 

so visual filter 31 to filter the previously conrputed output 
pixel values together with each of the possible output 
levels for the current pixel 32 to compute the visually 
perceived output value for each possible output level. 
For a binary output device there will be two possible out- 

55 put levels corresponding to a blad^ or a white pixel, but 
this method can be extended to multi-level output 
devices by considering more than two possible output 
lev^s. The output pixel value is chosen which gives the 
smallest difference between the continuous-tone input 
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value and the visually perceived output level. The result- 
ing error is then calculated for the current pixel by a dil- 
ference signal generator 33 which computes a 
cirfference signal between the continuous-tone input 
value and the visually perceived output value 34. As s 
with the conventional en-or diffusion technique, this dif- 
ference signal is then weighted by a series of error 
weights Wj j using a weighted error generator 35, and is 
added 36 to the continuous-tone input values of nearby 
pixels which have yet to be processed to form modified w 
continuous-tone input values^ The causal visual filter 
function is computed from the frequency response of 
the human viajal system as described by Sullivan. An 
example of a 4x7 causal visual f flter function is shown in 
FIG. 5. The array element 40 is used to weight the pos- 75 
sible output levels for the present pixel, and the remain- 
ing array elements are used to weight the nearby output 
values which have been previously computed. (This 
causal visual filter function is a normalized version of 
the one disclosed by Sullivan.) A similar 8x15 causal 20 
visual filter function is shown in FIG. 6. The size of the 
causal visual filter function array is related to the param- 
eters such as viewing distance ard sample spacing on 
the document in the equations given by Sullivan. In this 
case, the array element 42 is used to weight the possi- 25 
bie output levels for the present pixel, and the remaining 
array elements are used to weight the nearby output 
values which have been previously computed. FIG. 7 
shows a sample image generated using this technique 
with the 8x1 5 causal visual f flter function. The worm arti- 30 
facts 18A-C which were shown in FIG. 3 have been 
eliminated. 

Although visual error diffusion has the result of 
reducing the appearance of worms to a large degree 
without the introduction of undesirable noise, a side 35 
effect of this method is that artifacts are generated near 
edges, lines, and isolated pixels in the image. One type 
of artifact can be refen-ed to as "ghost pixels." Examples 
of ghost pixels artifacts 20A-C can be observed around 
the isolated pixels and the edges in the Image shown in 40 
FIG. 7, The origin of these gho^ pixel artifacts can be 
traced to the fact that the selection process is compar- 
ing the visually perceived output image to the unfittered 
input image. So, for example, consider the case where 
the input image contains an isolated black pixel on a 45 
white background. The visual error diffusion technique 
will place white pixels in most of the area conresponding 
to the white Ijackground. but when it comes time to 
make a selection for the black pixel it will visually filter 
the output pixels in the surrounding area. Since most of so 
these pixels will be white, the visually perceived output 
level for a white output pixel will be white, but the visu- 
ally perceived output level for a black output pixel will be 
light gray rather than blade Since the light gray visually 
perceived output level will be closer to the desir^ black ss 
input level than the white visually perceived output level 
, the correct black output level will be chosen. However, 
the resulting error signal corresponding to the difference 
between the black input level and the light gray visually 



perceived output level will be quite large. This error will 
then be propagated to the surrounding continuous-tone 
input pixels. As a result when the nearby white back- 
ground pixels are processed, one or more of these pix- 
els may be rendered as black output pixels. 

The artifacts can be even more severe on edges 
between uniform patches of different gray levels where 
anisotropic overshoots are formed. Although this artifact 
has a sharpening effect on the image, it is not necessar- 
ily desirable because the degree of sharpening can not 
be controlled independent of the worm reduction, and 
additionally the effect is anisotropic so that the apparent 
sharpening is not symmetric. Examples of the aniso- 
tropic sharpening of edges 19A and B can be seen In 
FIG. 7. As with the ghost pixel artifacts, the origin of tfie 
anisotropic sharpening artifacts can be traced to the fact 
that the selection process is comparing the visually per- 
ceived output image to the unfiltered input image. 

Ideally, it is desirable to conrtoine the performance 
of the error diffusion algorithm of FIG. 1 in areas of the 
image that contain high spatial frequency content 
(edges), and the performance of the visual enror diffu- 
sion algorithm in areas of the image that contain 
smoothly varying, low spatial frequency content 

Summary of the Invention 

It is an object of the present invention to eliminate 
the above-discussed artifacts while still providing worm 
reduction. 

This object is achieved in an error diffusion method 
suitable for producing an output image from an input 
image having a set of digitized continuous-tone pixels, 
comprising the steps of: 

a) corrputing an image activity signal; 

b) computing a set of activity weights from the 
image activity signal; 

c) computing a filtered input value for a digitized 
continuous-tone input pixel responsive to the activ- 
ity weights; 

d) computing a filtered output value responsive to 
ttie activity weights for each of the possible output 
levels; 

e) sdecting the output level in response to the fil- 
tered input value and the filtered output value for 
each of tiie possible output levels according to an 
error criterion; 

0 determining an error signal between the filtered 
input value arrd the filtered output value for the 
selected output level; and 

g) weighting tiie error signal and adjisting the fil- 
tered input values for nearby pixels which have not 
yet been processed. 

Advantaoes 

The present invention has the advantage that it 
minimizes the appearance of worm artifacts in the out- 
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put image whfle eliminating the edge artifects associ- 
ated with the prior art visual error diffusion algorithm 
disclosed by Sullivan. 

The present invention also has the advantage that it 
adapts to the local scene content by selecting a set of s 
operating parameters that produces the best image 
quality for that particular type of scene content. 

Another advantage of the present invention is that it 
pQ-mits the amount of sharpening to be adjusted inde- 
pendent of tile amount of visual blur. Acklitionally. the io 
sharpening effect can be made to be isotropic. 



FIG. 20 is a logic diagram showing a mettoi 
according to the present invention utilizing a single 
error buffer; 

FIG. 21 is a logic diagram showing a method 
according to the present invention applied to each 
channel of a color image; and 
FIG. 22 is another logic diagram showing another 
method according to the present invention appliej 
to each channel of a color image. 

Detailed Description of the Invention 



Brief Descrbtion of the Drawings 

FIG. 1 is a logic diagram of a prior error diffusion is 
process; 

FIG. 2 shows typical sets of error weights L^ed in 
the FIG. 1 arrangement; 

FIG. 3 depicts a sample image generated in 
accordance with the method of FIG. 1 which con- 20 
tains worm artifacts; 

FIG. 4 is another logic diagram of a prior error diffu- 
sion process; 

FIG. 5 depicts a causal visual filter array: 
FIG. 6 depicts another causal visual filter array; 25 
FIG. 7 depicts a sample image generated in 
accordance with the method of FIG. 4 which con- 
tains anisotropic edge sharpening and ghost pixel 
artifacts; 

FIG. 8 is anotiier logic diagram of an error diffusion 30 
process; 

FIG. 9 depicts a sample sharpening filter which can 
- be used in ttie logic of FIG. 8; 

FIG. 10 depicts a sample image generated in 

accordance with the method of FIG. 8 which elimi- 35 

nates the artifacts of FIGS. 3 and 7, but contains 

artifacts in high contrast image regions; 

FIG. 11 is a logic diagram showing a method 

according to the present invention; 

FIG. 1 2 is a logic diagram showing an input filtering 40 

operation according to the preserrt invention; 

FIG. 13 Is a sample graph describing a function for 

generating a waghting vector Q according to the 

present invention; 

FIG. 14 is another sample graph describing a tunc- 45 
tion for generating a weighting vector Q according 
to the present invention; 

FIG. 15 is a logic diagram showing a weighting of 
the error buffers according to the present irtvention; 
FIG. 16 is a logic diagram showing an output filter- so 
ing operation according to the present invention; 
FIG. 1 7 is a logic diagram showing another output 
filtering operation according to the present inven- 
tion; 

FIG. 1 8 is a logic diagram showing distributing error 55 
into separate error buffers according to the present 
invention; 

FIG. 19 depicts a sample image generated in 
accordance with the present invention; 



As set forth in the above referenced copending 
application, a modification to the visual error diffusion 
algorithm that corrects for these undesirable artifacts 
has been disclosed in which a visual filter function is 
applied to both the input and output images during the 
process Of determining the output value for the current 
pixel. The error that is propagated to nearby continuous- 
tone input pixels that have not been processed yet is 
calculated by the difference between the visually per- 
ceived input value and the visually perceived output 
value. This invention has tiie advantage that it mini- 
mizes the appearance of worm artifacts in the output 
image while eliminating the edge artifacts associated 
with ttie prior art visual error diffusion algorithm. A block 
diagram of the improved visual error diffusion algorithm 
is shown in FIG. 8. and a sample output image is shown 
in FIG. 1 0. One drawback of the FIG. 8 metiiod is that in 
order to restore the desirable appearance of sharpness 
to the image, a sharpening pre-filter 68 is required. An 
example sharpening pre-filter. shown in FIG. 9. was 
used in creating the image of FIG. 10. Thus, the edges 
in FIG. 10 appear as sharp as the edges in FIG. 7, 
which was processed with the visual error diffusion 
method of Sullivan, but witii the ackJed advantage of 
being isotropic in nature. While this improvement pro- 
duces a higher quality output than the previous meth- 
octe. ft does introduce undesirable noise 21A-B into 
certain areas of the image that are high in contrast and 
spatial frequency, such as the vertical and horizontal 
bars near the center right of the sample image of FIG. 
10. 

In the present invention, shown schematically in 
FIG. 1 1, a pixel located in row i, column j of the input 
image is denoted by y(i,D. An activity detector 80 uses 
naghboring pixels in the vicinity of y(iii) to determine an 
activity signal. TTie activity detector 80 may take many 
forms, including a local range detector, local variance 
estimator, convolution filter, or edge detection operator 
such as a Sobel or Prewitt edge detection filter. In one 
version, the activity detector 80 uses a local range 
detector to compute the activity signal as the difference 
between the maximum and minimum input pixel values 
of neighboring pixels in the vicinity of y(ij). The neigh- 
boring pixels may include the current pixel, adjacent pix- 
els, and/or other nearby pixels. Another activity detector 
80. in accordance with the inver^on, uses a local vari- 
ance estimator to compute the activity signal as the sta- 
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tistical variance of the input pixei values of neighboring 
pixelsin the vicinity of y(i.j)- Stfll another activity detector 
80 uses a convolution filter to compute the activity sig- 
nal as the convolution of the convolution filter and the 
input pixel values of neighboring pixels in the vicinity of 
y(i,j). An activity detector using an edge detection oper- 
ator computes the activity signal as the sum of convolu- 
tions of the input pixels of neighboring pixels in the 
vicinity of y(i,j) with the edge detection filters. 

The value of the activity signal is used as an index 
to an activity function 82 which computes a weighting 

vector Q = {qo. Qi Qk* -•*^n-i}- '"^ ^® prefen-ed 

embodiment of this invention, the activity function 82 is 
implemented as a lookup table tiiat is indexed by the 
activity signal. The input image is also processed 
through an input filter processor 81 , which digitally filters 
the input image in the vicinity of y(i J) using filtered signal 
generators 81 A and a bank of N separate filters FqO.O 
through FM.i(ij) as shown in FIG.12. Also in F1G.12. the 
weighting vector Q is used to weight the N filtered ver- 
sions of the input image using multipliers 81 B. which are 
then added together using an adder 81 C to produce a 
weighted filtered input signal 2(i.j). The filtered input val- 
ues only need to be computed for cases where the 
weights are non-zero. 

An example of an activity function 82 that can be 
used to switch between two input filters Fo(i.D and F^ilii 
is shown in FIG. 13, where |e| is the value of the activity 
signal. In this particular exarrple. the value of the activ- 
ity signal ranges from 0 to 255. For values of the activity 
signal between 0 and 10, Q = {0,1}. and ttie input filter 
F^(i,j) will be used. For values of tiie activity signal 
between 10 and 255. Q = {1 .0), and the input filter Fo{i,i) 
will be used. In this way. tine method of the present 
invention adapts the image processing t>ased on the 
value of the activity signal. This feature can be general- 
ized so that a different input filter is selected for each of 
N types of local scene content in the input image, as dis- 
criminated by the activity detector. Additionally, the use 
of tiie filters Fo(i,j) through F^^iO j) can be blended in 
regions of tine image where it is unclear exactly what 
type of scene content is present. An example of an 
activity function 82 which will accomplish tiiis is shown 
in FIG. 14. 

At iSne output of the input filter processor 81, it is 
desired that the average value of z(i.i) be equal to the 
average value of y(i,j) for uniform regions of the input 
image. One way of satisfying this average preserving 
relationship is to use filters Fo(i,j) through Ff^.ip.j) tiiat 
are normalized (i.e.. tiiat satisfy tiie relationship 
5:,2jF(i,j) = 1) . and impose the consti-airrt that Q satis- 
fies the relationship £ jQj = 1 . 

After the weighted filtered input signal z(i.j) is com- 
puted, tfie errors made by processing previous pixels 
are weighted by the en-or weights W(i,j) shown in FIG. 2 
and by the weighting vector Q using a weighted en^or 
generator 83 (FIG. 11). The details of the weighted error 
generator 83 are shown in FIG. 15. in which error values 
stored in N separate en-or buffers boO.i) through bN.i(i.j), 



corresponding to the N types of local scene content, are 
individually weighted by W(i.j) using filtered enor gener- 
ators 83A, then weighted by the the weighting vector Q 
t^ng multipliers 83B, and ackjed together using an 

5 acWer 83C to produce an error signal e(i,j). In this way. 
the errors generated in areas of the image that contain 
one type of scene corrterTt will not propogate to an area 
of the image that contains a different type of scene con- 
tent This type of s^Darate error buffering is desirable to 

10 avoid artifacts that can be generated when the image 
contains a sharp boundary between areas of different 
scene content, such as the boundary between the 
crosshair patches and the horizontal and vertical sti-ipe 
patdies near the cerrter right of FIG. 10. 

75 Returning to FIG. 1 1 , the error signal e(i,j) is added 
to the weighted input signal z(i,j) by an adder 88 to com- 
pute the desired signal d(i,j). An output filter processor 
89 conputes a set of M weigfited filtered output values 
Pq through p^.^. This is done by digitally filtering the 

20 previously computed output values o(i.j) together with 
each of the M possible output levels 85 for the current 
pixel (lo through If^i) using filtered signal generators 
89A and the bank of N filters Go(i,j) tiirough GN.i(i.j) as 
shown in FIG. 16 to compute filtered output values. For 

25 each of the M possible output levels, the filtered output 
values corresponding to the N filters are weighted by the 
weighting vector Q using multipliers 898, and added by 
an adder 89C to generate the weighted filtered output 
values, as shown in FIG. 16. The filtered output values 

30 only need to be computed for cases where the weights 
q^ are non-zero. 

Returning again to FIG. 11, once the weighted fil- 
tered output values po through p^-i have been com- 
puted, a selector 86 chooses the output level Iq through 

35 according to an error criterion. One such error crite- 
rion is to minimize the difference between the desired 
signal d(i.j) and tiie weighted filtered output values po 
ttirough Pm-i- 

Once a selection of an output level has been made. 

40 an output filter processor 87 digitally filters the output 
signal o(i,j) using filtered signal generators 87A and the 
bank of N filters Go(i.i) through GN.i{i,j), weights the 
results by the weighting vector Q using multipliers 878. 
and adds the weighted values together using an adder 

45 87C as shown in FIG. 17. An error signal generator 90 
(FIG. 11) then generates an error signal by taking the 
difference between the value generated by tiie output 
filter processor 87 and tiie desired signal d(i.j). The error 
signal is then distributed into N enor buffers bo through 

50 bM.i by an error distributor 84 which weights the error 
signal by the weighting vector Q using multipliers 84A 
as shown in FIG. 18. The values stored in the N error 
buffers bo through b|sj.i are then weighted by the 
weighted error generator 83 and added using the adder 

55 88 to weighted filtered input pixels which have not yet 
been processed, as described earlier. 

A sharpening pre-filter can be applied to the input 
image in order to restore the desirable edge sharpness. 
If the sharpening is implemented by a device which per- 
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forms a convolution, then it is possible to incorporate the 
sharpening pre-f ilter directly into the input filtering oper- 
ation 81 by convolving the sharpening pre4ilter with the 
filters Fo(i.j) through FM.i(ij) to arrive at a new set of 
input fOters. It may be desirable to only apply the sharp- 5 
erring filter to portions of the innage which contain large 
amounts of activity. This can be accomplished by only 
incorporating the sharpening pre-filter into a subset of 
the filters Fo(i.i) through F^.^pj) which correspond to 
large values of the activity signal. 70 

A sample image generated using the method of the 
present invention is shown in FIG. 19. The activity func- 
tion shown in FIG. 13 was use6 to switch between 
image regions having high and low activity signal val- 
ues. The activity detector 80 used in this example was a 15 
local range detector In regions of low activity, the input 
filter Fi(i.j) used was a delta function 6(i.j), where 6(i.j) is 
defined as 1 for i=j=0 and 0 otherwise. As a result, the 
input filter operation reverts to a null operation and z(i,j} 
is equal to y(l.j). The output filter used for low activity 20 
regions G^O.j) was set equaJ to the causal visual filter 
V(ij) shown in FIG. 6. 

In the high activity regions, the input filter Fo(i»i) 
used was the sharpening fOter 68 shown in FIG. 9. The 
output filter GoOJ) used was the delta function 5(iJ). 2s 

For values of the activity signal between 0 artd 10. 
? Q = {0,1}. and the visual error diffusion method of FIG. 
4 is used. Alternately, for values of the activity signal 
between 10 and 255, Q = {1.0}. and the conventional 
error diffusion of FIG. 1 with a pre-sharpening step is 30 
used. In this way, the method of the present invention 
adapts the processing based on the value of the activity 
signal. 

Another embodiment of the invention enrtploys a 
single buffer to store the error signals generated by the 35 
error signal generator 90, as shown in FIG. 20. Compar- 
ing this with the metiiod shown in FIG. 11 , it can be 
seen that the N separate error buffers bo through b|v|.^ 
have been eliminated arKl the error signal feeds directiy 
into a weighted error generator 91. The weighted error 40 
generator 91 then computes a weighted error signal 
e{i.j) using the weighting vector Q and the error signal 
from the error signal generator 90. 

For color images, the method described above can 
be applied to each channel of the innage as shown in 45 
FIG. 21 . For example, if an image consists of red. green, 
and blue pixels, the improved adaptive error diffusion 
' method can be applied via block 70 to the red pixels, 
and again to the green pixels and the blue pixels by 
block 71 and 72, respectively. Similarly, the improved so 
adaptive error diffusion method can be applied to each 
channel of a 4-color cyan, magenta, yellow and black 
image. It may be desirable to use different input filters, 
output filters, activity lookup tables, error weights, activ- 
ity functions, and/or sharpening filters can be used for ss 
each color channel. 

Another method of processing color images 
according to the present invention is to process all of the 
color channels simultaneously, as shown in FIG. 22. It 



can be seen that the logic of FIG. 22 is the same as the 
logic of FIG. 1 1 , with the exception that the signals that 
are communicated between logical blocks are now vec- 
tor signals, as opposed to scalar signals discussed in 
ttie method of FIG. 11. In the metiiod of FIG. 22. the 
input continuous-tone pixel values for row i column j of 
the color channels are combined into a color vector 
Y(i,j). A vector activity detector 100 determines an activ- 
ity signal from Y(i.j). and the activity function 82 deter- 
mines a weighting vector Q from the activity signal- The 
input color vector Y(i,j) is then filtered by a vector input 
filter processor 101. which generates a filtered input 
color vector Z(i,j). A weighted color vector error E{i,j) is 
computed by a weighted color vector error generator 
102 and added to Z(i,j) by a vector adder 108. resulting 
in a desired color vector D(i,j). A vector output filter proc- 
essor 109 determines a set of weighted filtered output 
color vectors Pq through P^.^ by filtering previous out- 
put color vectors O(i.j) along with possible output color 
vectors 105 (Lq through L^.^) using a set of output filters 
Gq through G^.i and the weighting vector Q. A selector 
106 chooses an output color vector according to an 
error criterion. One such error criterion is the minimum 
vector distance between the desired color vector D(l,j) 
and the weighted filtered output color vectors Pq 
through P^.^. Once an output color vector has been 
chosen, a vector output filter processor 107 filters tine 
output color vector 0(i,j) using a set of output filters Gq 
tiirough G^^.^ and the weighting vector Q. A cojor vector 
enor signal generator 1 10 computes a color vector error 
signal as the vector difference between the output of the 
output filter processor 107 and the desired color vector 
D{i,j). The color vector error signal is then distributed 
into N vector error buffers Bq through B^.^ using a vec- 
tor error distributor 1 03 which is responsive to the 
weighting vector Q. The N vector error buffers Bq 
tiirough 8^.^ are then used by the weighted color vector 
enror generator 1 02 to compute a weighted color vector 
error E(i.j) for use in processing future color vectors. 

The invention has been described in detail with par- 
ticular reference to certain preferred entDOdiments 
thereof, but it will be understood tiiat variations and 
modifications can be effected within the spirit and scope 
of the invention. 
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8. The method of claim 1 wherein the weightir^ the 
error signal step includes weighting the error signal 
as a function of the activity weghts. storing such 
weighted error signals, and adjusting the stored 
weighted error stgral as a function of a set of error 5 
weights and the activity weights. 

9. The method of claim 1 wherein the number of activ- 
ity weights is 2. 

10 

10. The method of claim 13 wherein the input filter cor- 
responding to the first activity weight is a delta func- 
tion, and the input filter corresporKiing to the 
second activity weight is a delta function. 
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(54) Adaptive error diffusion method 

(57) An error diffusion method suitable for produc- 
ing an output image from an input image having a set of 
digitized continuous-tone pixels is disclosed. The 
method includes the steps of computing an image activ- 
ity signal; computing a set of activity weights from the 
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ble output levels. The method further includes selecting 
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put levels according to an error criterion; determining an 
error signal between the filtered input value and the fil- 
tered output value for the selected output level; and 
weighting the error signal and adjusting the filtered input 
values for nearby pixels which have not yet been proc- 
essed. 
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